WELL-POSEDNESS RESULTS FOR TRIPLY NONLINEAR 
DEGENERATE PARABOLIC EQUATIONS 

B. ANDREIANOV, M. BENDAHMANE, K. H. KARLSEN, AND S. OUARO 

Abstract. Wc study the wcU-poscdncss of triply nonlinear degenerate elliptic- 
parabolic-hyperbolic problem 

b{u)t — div a{u, V4){u)) + ip{u) = /, u\t=o = "0 

in a bounded domain with homogeneous Dirichlct boundary conditions. The 
nonlinearities b, </> and ip are supposed to be continuous non-dccroasing, and the 
nonlinearity a falls within the Leray-Lions framework. Some restrictions are 
imposed on the dependence of a{u, V(^{ji)) on u and also on the set where ip de- 
generates. A model case is o(u, V0(u)) = f(6(M), i/'(m), </)(n))-|-fc(M)ao( V</)(m)), 
with 4> which is strictly increasing except on a locally finite number of seg- 
ments, and Oq which is of the Leray-Lions kind. We are interested in existence, 
uniqueness and stability of entropy solutions. If b = Id, we obtain a general 
continuous dependence result on data mq,/ and nonlinearities b,ip,<j>, a. Simi- 
lar result is shown for the degenerate elliptic problem which corresponds to the 
case of 6 = and general non-decreasing surjective tji. Existence, uniqueness 
and continuous dependence on data uq , / are shown when [b + ■(/)] (R) = R and 
(f) o [b + ij}]^^ is continuous. 
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1. Introduction 

1.1. Problem and assumptions. In this paper we consider problems under the 
general form 

( dth{u) + div f(6(w), ^(w), (j){u)) - div a(u, V<?!)(u)) + i'{u) = /, 
(P) \ in Qt = (0, T) X 

[ u|t=o = wo in fi, u = on(0,T)x9ri, 

where u : (t, x) G Qt — *■ R. is the unknown function, T > is a fixed time, 
f2 C is a bounded domain with Lipschitz boundary 9f2. 
We assume 

the functions &, '0, ^ : R i— > R are continuous nondecreasing, 
normalized by the value zero at the point zero. 

We require the following technical assumption on </>: 



{Hi) 



there exists a closed set £' C R such that cf) is strictly increasing on R \ i?, 
and the Lebesgue measure meas </>(£') is zero; 



and, moreover, 

(i?3) liminfexo meas {G^)/e < +oo, where G"^ := {z e R | dist (z, (f>{E)) < e}. 



Notice that since is continuous and strictly monotone on R \ P, the set G := 
(j){E) is also closed. 

Remark 1.1. (i) Hypotheses {H2),{H^) are trivially satisfied if (p is a strictly in- 
creasing function. In case (j) has a finite number of segments on which it keeps 
constant values, E is just the union of all these "flatness segments", and {H2),{H3) 
are satisfied. 

(ii) Property {H2) is still true if (j) is locally absolutely continuous. In general, 
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the set of discontinuity points of (f>^^ is not closed, and its closure can be large 
(this is the case, e.g., if (j) is the "Cantor stairs function"). Thus {H2) is a re- 
striction, although it is fulfilled in most of the practical cases. Property is a 
further restriction. Indeed, consider the following example. It is easy to construct 
a Lipschitz continuous non- decreasing function (j) such that G ~ 4>{E) is equal to 
{0} U G N}. A straightforward calculation shows that for e = \/n, 
contains the whole interval [0, l/(2n^/'^)]; in this case meas(G"^)/e is of order e~^^^ 
and gets unbounded as £ — *■ 0. 

The initial function uq : 51 ^ R and the source / : Q ^ R are assumed to fulfill 
uo G L'^{Vl); f is measurable such that 



fit, 



eL°°(r!)fora.e. tG(0,T)and / •)IIl-(o) < 

"'0 



Furthermore, the following condition (automatically satisfied in the case 6(R) = R) 
is needed : 



in the case 6(+oo) 7^ +00 (rcsp., 6(— 00) 7^ ~oo) one has 
^{+00) = +00 and /+ e L'^iQr) 

( resp., 00) = —00 and € L°°{Qt) )■ 



Assumptions {H4),{Hz) are imposed to limit our study to bounded solutions of (P). 

Note that in view of {Hi) and (i/s), we are assuming at least that (6+?/')(R) = R. 
An important part of the paper in devoted to the case 6(R) = R. If, in addition, 
6 is bijective, then performing a change of the unknown one puts the problem into 
the doubly nonlinear framework with b = Id. 

Our continuous dependence result for problem (P) (in which we perturb both 
the data and the nonlinearities) concerns the case where the structure condition 

[Hstr) b{r) = h{s) ^ 4>{r) - 4>{s) 

is satisfied. This result implies the existence of solutions for (P), by reduction 
to non-degenerate problems. Assumption (Hstr) is trivially satisfied in the case 
b = Id. If {Hstr) fails, the convergence of approximate solutions to (P) is known 
only for a particular monotone approximation method developped by Ammar and 
Wittbold [4]. This approach leads to an existence result which bypasses {Hgtr)', 
but interesting issues (such as proving convergence of numerical methods for (P) 
without requiring the structure condition {Hstr)) remain open. See Benilan and 
Wittbold [J 4] for a thoroughful discussion of the role of the structure condition for 
a simple model one-dimensional case dtb{u) -I- {\{u))x = u^x- 
Furthermore, 

(-ffg) the function f : R x R x R ^ R^ is assumed merely continuous. 

Notice that under the structure condition {Hstr), the dependency of f on (f>{u) can 
de dropped. Whenever it is convenient (and in particular, in the case where b is 
bijective), we use the notation f(-) := f(6(-), ■(/;(•), 0(-)). 

The function a : R x R^ R^ is assumed to satisfy the following conditions : 

{Hj) a is continuous" on R x R^, and a{r, 0) = 0; 



^In the sequel, we will abusively denote the latter quantity by 



llil{0,T;L°°{n))- 
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iH9) 



a(r, •) is monotone, i.e. 

(a(r, - a(r, 77)) • - r?) > for all ^, 77 G and all r G R; 

a(r, •) is coercive at zero, i.e., there exist p € (1, +00) and C G C(R; R+) 
such that a(r, ' C > c(Fyl?l^ f G R^ and aU r G R; 

the growth of a(r, ^) is not greater than i-C, 
there exists C G C(R; R+) such that 

|a(r,OI < (^(0(1 + \^\P~^) for all ^ G R^ and aU r G R. 

It follows from {Hj)-{Hio) that for all r, the operator w — div a(r, Vw) is an 
operator acting from WQ^'^(r2) to W^^'^ (fl), where p' = Since the work of 

Leray and Lions [29] , this assumption became classical. It can be generalized to the 
framework of Orlicz spaces (see the works of Kacur [27] and those of Benkirane and 
collaborators (sec e.g., [16], [15]), and even to more general coercivity and growth 
assumptions. We refer to Bendahmane and Karlsen [10, 11] for the case of the 
anisotropic p-Laplacian. In the case of dimension = 1, very general coercivity 
assumption lim^^-too — (uniformly for r bounded) can be considered. 

In this framework, the well-posedness for (P) was already established by Ouaro and 
Toure [37, 38, 39] and Ouaro [30] (see also Benilan and Toure [13]); notice that some 
essential arguments of these works are specific to the case TV = 1. Notwithstanding 
the above generalizations, the classical Leray-Lions assumptions arc sufficient for 
us to illustrate the arguments of the existence proof for (P). 

The relevant technical assumption in order to have uniqueness is 



(^11) 



there exists C G C(R^;R+) such that 

(a(r, - a{s, r,)) ■ (C - r/) + C(r, s){l + \^\p + MniHr) - -^(s)! > 
for all ^, ?7 G R^ and all r, s G R such that 
the segment which lies between r and s 
does not intersect the exceptional set E. 



This assumption goes along the lines of Carillo and Wittbold [ ' i ] and combines 
the monotonicity condition (iJg) with a kind of Lipschitz continuity assumption on 
a(-, ^) o on the connected components of R \ i?. 

Remark 1.2. Notice that E is a set of the values of u that lead to a{u, V0(w)) 
being zero. Indeed, since meas 4>{E) = 0, we have V0(w) =0 a.e. on the set where 
u G E; then we have a(u, V(/)(u)) = a(ii,0) = 0, regardless of the exact value of 
ue E. 

Remark 1.3. Notice that we do not assume the structure condition 

Hr)=c^{s) => a(r,C) =a(s,e) G R^. 

This means that a(-,^) can be discontinuous with respect to (j){r); the set of discon- 
tinuities is included in 4>{E) which, by {H2), is a closed set of measure zero. This 
technical assumption is needed to be able to "cut off" the discontinuity set. 

One can also consider a(r, ^) which is discontinuous in r. E.g., take the case of 
'^{''^tO — '^iikir),S.) with fc(-) piecewise continuous. Thanks to Remark 1.2, it is 
reduced to our setting by a change of unknown function u into v such that u ~ p{v) 



^ the assumption of continuity of o in the first variable can be relaxed: see Remark 1.3 
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with p non-strictly increasing, chosen so that k{p[z)) = k{z) with k(-) continuous. 
Indeed, such change of the unknown preserves assumption (Hgtr)- 

Let us mention that the assumptions {Hi),{H5) and (Hg), (-ffio) can also be 
generahzed. Within the framework of "variational" solutions, one usually works 
with "bounded energy initial data" , i.e., with uq measurable and such that B{uq) < 
+00, where 

(1) B{z) := r <j,{s)dh{s) 

Jo 

is the function depending on h and (j) which we introduce following Alt and Luckhaus 
[1], and with relaxed growth and coercivity assumptions allowing for additional 
terms which depend on B{r) + |^(r)p + rip{r), these terms being controlled by 
means of a priori estimates (see, e.g., [1, 27, 5]). 

A more general framework is provided by the one of renormalized solutions (see 
[18, 19, 21, 4, 17, 10] and the references therein). Nonetheless, the assumptions 
{H4), {H5), (Hg), and (J^io) are sufficiently weak to provide the starting point for 
the well-posedness theory for renormalized solutions of (P). Indeed, the uniqueness 
proof for renormalized solutions of (P) remains essentially the same as the one of 
Carrillo and Wittbold in [21], and the existence result is most easily obtained 
using bi-monotone sequences of bounded entropy solutions, following Ammar and 
Wittbold [4] and Ammar and Redwane [3]. 

1.2. The notion of solution and known results. Problem (P) is of mixed 
elliptic-parabolic-hyperbolic type, and thus combines the difficulties related to non- 
linear conservation laws with those related to nonlinear degenerate diffusion equa- 
tions. We refer to Kruzhkov [28] and to Leray and Lions [29], Lions [30], Alt and 
Luckhaus [1], Otto [35] for the fundamental works on these classes of equations, 
respectively. 

One consequence is that the notion of weak solution (sometimes called "varia- 
tional solution") generally leads to non-uniqueness, unless </)(•) is strictly increasing. 
The notion of entropy solution we use is adapted form the founding paper of Carrillo 
[20] , which extends the classical framework of entropy solutions to scalar conserva- 
tion laws to the case of problem (P) with the linear diffusion a{u, V0(u)) = Vc/)(u). 
The uniqueness arguments of [20] were adapted by Carrillo and Wittbold [21] to 
the case of a nonlinear Leray-Lions diffusion operator of the form a{u, Vu) corre- 
sponding to (j) — Id. The case of 6 = Id and of a diffusion operator of the form 
a{(j){u), V0(u)) is similar; one particular case is considered in Andreianov, Ben- 
dahmane and Karlsen [(>]. Both frameworks are sometimes referred to as "doubly 
nonlinear" . Notice that a new definition of an entropy solution, suitable also for 
the case of doubly nonlinear anisotropic diffusion operators, was used in a series of 
works by Bendahmane and Karlsen (see [10, 11] and references therein); the issue of 
existence in this general anisotropic framework is still open. The case of triply non- 
linear problems of the form (P) has been first addressed by Ouaro and Toure (see 
[39] and the references therein) and Ouaro [3G] . Well-posedness results are obtained 
in dimension = 1, under very general coercivity conditions; see also the works of 
Bcnilan and Toure ([13] and the references therein). The multi-dimensional elliptic 
analogue of (P) was recently addressed in the work Ammar and Redwane [3], in 
the framework of renormalized solutions; their approach is quite similar to the ours. 
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but the proofs of [■'>] require a special structure of the nonlinearity 0. We bypass 
the difficulties of [3] using two observations in Section 3 (see Lemmas 3.2,3.3). 

In most of the works cited hcreabove, the homogeneous Dirichlet boundary con- 
ditions were chosen. One should bear in mind that, unless (f) is strictly increasing, 
the boundary condition is also understood in an entropy sense (see e.g. Bardos, 
LeRoux and Ncdelec [0], Otto [14], Carrillo [2(J]). Focusing on the homogeneous 
boundary condition is a simplification which seems to be not merely technical. A 
partial extension of the Carrillo's techniques to inhomogeneous boundary data can 
be found in Ammar, Carrillo and Wittbold [2]. Different techniques for the inho- 
mogeneous problem, based on the weak trace framework introduced by Otto [34] 
and developed by Chen and Frid [22], were used by Mascia, Porretta and Terracina 
[32] and Michel and Vovelle [33]. A related, though somewhat more straightforward 
approach was attempted by Andreianov and Igbida [7]. Notice that in all cases, a 
technique of "going up to the boundary" is preceded by obtaining the fundamental 
weak formulations and entropy inequalities "inside the domain" . In the present 
paper, we also make the simplest choice of the homogeneous Dirichlet data, and 
focus on deriving the entropy inequalities. 



1.3. Main techniques and the outline of the paper. Our main concern is 
the existence for (P) (and more generally, the continuous dependence result with 
respect to perturbations of the data and the nonlinearities) in the case b ^ Id 
(or, more generally, under the structure condition {Hgtr))- We extend the ar- 
guments of Andreianov, Bendahmane and Karlsen [6] developed for the case of 
a{u, V0(m)) = a{\/(j){u)); the main difficulty stems from the fact that we do not 
assume the structure condition of Remark 1.3, so that a(0^^(-),^) can be discon- 
tinuous. 

In order to use the weak convergence while passing to the limit in the nonlinear 
diffusion term in (P) , we use a version of the classical Minty-Browder monotonicity 
argument. We use cut-off functions to focalize on the intervals of the strict mono- 
tonicity of (j). We then show that the complementary of this set can be neglected, 
thanks to a particular a priori estimate with the cut-off function which localizes on 
a neigbourhood of the "exceptional set" E of the values of u introduced in {H2) 
(see Remark 1.2 and Lemma 3.2). 

In order to deal with the convection term in (P), we use the technique of non- 
linear weak-* convergence to a measure-valued solution (as considered by Tartar, 
DiPerna, Szepessy, Panov), or more exactly, to an entropy process solution as de- 
veloped by Gallouet and collaborators (see [25, 23, 24] and references therein). 
Considering entropy process solutions is a purely technical issue, since the unique- 
ness proof also contains their identification to entropy solutions. 

The chain rule arguments of Lemmas 3.3,3.4 permit to separate the two afore- 
mentioned weak convergence arguments, the one for the diffusion term and the one 
for the convection term. 

The uniqueness of an entropy solution is shown under the assumption (Hu), with 
the help of {H3) and the estimate of Lemma 3.2; we follow Carrillo [20], Carrillo 
and Wittbold [21], Eymard, Gallouet, Herbin and Michel [24] and Andreianov, 
Bendahmane and Karlsen [G]. 
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While the full continuous dependence result strongly relies upon the structure 
condition {Hstr), the stability of entropy solutions to (P) with respect to the per- 
turbation of the data (uq, /) is shown under the weaker structure assumption 

{HUr) {b + i'Kr) = (& + ^ 0(r) = 

Let us stress the fact that our argument using {H'^^^) in not "robust", in the sense 
that it cannot be directly adapted to the proof of convergence of various kinds of 
approximate solutions. In order to address the question of convergence of numer- 
ical approximations of (P), the stronger structure assumption (Hgtr) still seems 
essential. For finite volume approximations of the doubly nonlinear problem (P) 
with 6 = Id, a convergence proof using (Hstr) is given in [(>]. 

Let us give an outline of the paper. We start with definitions and the formulation 
of the main results in Section 2. In Section 3, we give the key ingredients of 
our techniques. Section 4 concerns the adaptation of the standard uniqueness, 
contraction and comparison result for entropy and entropy- process solutions of (P). 
Section 5 contains the a priori estimates for solutions. In Section 6, we assume 
that b is bijective (or, more generally, conditions {Hstr),{H5) are satisfied); we deal 
with the convergence of solutions to problems {Pn) with perturbed coefficients. 
Section 7 is devoted to the proof of the well-posedness result for (P). In Section 8 
we give existence, uniqueness and continuous dependence results for the related 
elliptic equation -I- div f(?/'(u), 0(u)) — divo(w, V0(w)) = sg L°°(r2). 

2. Entropy solutions and well-posedness results 

2.1. Entropies and related notation. As it was explained in the introduction, 
we need the notion of weak solution for (P) with additional "entropy" conditions. 
In order to use entropy conditions in the interior of Qt and, moreover, take into 
account the homogeneous Dirichlet boundary condition, following Carrillo [20] we 
will work with the so-called "semi-Kruzhkov" entropy-entropy flux pairs (?7^,ci^) 
for each c S R: 

Vt (z) = [z ~ c)+ , Tj- (z) = {z - c)- , 
q+(z) = sign+ (z - c) (f(z) - f(c)), q-(z) = sign- (z - c) (f(z) - f(c)). 
By convention, we assign (?7^)'(c) to be zero. Here (z— c)^ stand for the nonnegative 
quantities such that 2 — c = (z — c)+ — (z — c)~, but we denote 

z > c 



sign'*' 



(z-c) = (,y+)'(z) = { J; ^ < ^ sign- (z-c) = (,?-)'(z) = I 



— i, z < c 



At certain points, we will also need smooth regularizations of the semi-Kruzhkov 
entropy-entropy flux pairs; it is sufficient to consider regular "boundary" entropy 
pairs {'r]t,e,c{t,e) (^f- Otto [34] and the book [31]), which are W^'°° pairs with the 
same support as {rj^, q^), converging pointwise to (77^, q^) as e ^ 0. In particular, 
the functions 

signj'(z) = -min{z'^,e}, sign~(z) = -max{— z~,— e} 
e e 

will be used to approximate sign^(-) = (f?^ )'(•)■ 

Definition 2.1. For a function ip which is monotone on R, for all locally bounded 
piecewise continuous function 9 onR we can define (using, e.g., the Stiltjes integral) 



(/Sg : z G R i-^ / 6{s) dip{s) 
Jo 
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Moreover (see Lemma 3.1 below), there exists a continuous function (fe such that 
(feiu) = ip0{ip{u)). 

In the sequel, we denote by bf{-) the function z i-^ (vtYi^) db{s). 

Jo 

2.2. Entropy and entropy process solutions. For the sake of simplicity, we will 
in this paper work with bounded entropy solutions, i.e., we require that u S L°°{Qt) 
and put corresponding hypotheses on the data uo, f and functions b, tJj. Note that 
the boundedness assumption can be bypassed in the framework of renormalized 
solutions (see in particular Ammar and Redwane [■>] ) ; or in the framework of varia- 
tional solutions in the spirit of Alt and Luckhaus [1] (in this case, one has to replace 
the functions C(-) in assumptions (Hg), (i?io), (Hu) with a constant C; further 
changes are indicated in Remark 6.1). 

Definition 2.2 (entropy solution). An entropy solution of (P) is a measurable 
function u : Qt R satisfying the following conditions: 

(D.l) (regularity) u G L°°{Qt) and w = 0(u) G LP{0,T; W^'^in)). 
(D.2) For all ( G X>([0, T) x fl), 

b{u)dti + Ku) ■ - a(u, Vw) • - V(w)C ) dx dt 

(2) 

&(wo)e(0, ■)dx + I I dxdt = 



(D.3) For all (c,^ G x V{[0,T) xfl), ^ > 0, and also for all (c,C) G R x 
2?([0,T)xr!), e>0, 

bt{u)dt^ + qt{u) ■ \/^-{rjty{u)a{u,^w) ■ - {vtYi^W^)^^ dxdt 
bt{uo)m •) dx+ f f (vfYiu) f^dxdt > 0. 

n J JQt 

Remark 2.1. If in the above definition, u satisfies (D.l), if (2) is replaced by 
the inequality ">" (resp., with the inequality "<"), and if (D.3) holds with the 
entropies ry+ for c G R"*" (resp., with the entropies r/^ for c G R~J, then u is called 
entropy subsolution (resp., entropy supers olution) of{P). 

Remark 2.2. Following Alt and Luckhaus [\], we can rewrite the weak formulation 
(2) of (P) as follows: 

- b(u)\t=o — b{uo) and the distributional derivative dtb{u) satisfies 

dtb{u) G LP\0,T; W-^'P'm + L\Qt) 

in the sense 

'^{dtb{u) , = - / / Ku) dtC - / 6(^io)C(0, •) 

J JQt JU 

for all C G LP(0, T, W^'^in)) n L°^{Qt) such that dtC G L^{Qt) and C{T, ■) = 0; 

- equation (P) is satisfied in LP {0,T;W-^'P' (n)) + L^{Qt) ■ 

We denote by {■ , • ) the duality pairing between LP (0, T; W~^'P (fi)) + L^{Qt) and 
LP{0,T,W^'P{n))nL^{QT). 
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For technical reasons, it is convenient to introduce the notion of entropy pro- 
cess solution adapted from Eymard, Gallouet and Herbin [2:!], Gallouct and Hu- 
bert [25] and Eymard, Gallouct, Herbin and Michel [24]. This definition is based 
upon the so-called "nonlinear weak-* convergence" property (see e.g.. Ball [S] and 
Hungerbiihler [26]): 



(3) 



each sequence (uri)n of measurable functions 
admits a subsequence such that for all F E C(R, R), 



/ F{^{-,a)) da weakly in L^{Qt) 
Jo 



whenever the set {F{un))n is weakly relatively compact in L^(Qt), 



where the function jj, £ L°°{Qt x (0,1)) is referred to as the "process function". 
Notice that in the above statement, one also concludes that F{^{-, a)) independent 
of a whenever F(un) converges to F{u) in measure (see Hungerbiihler [2()]). 

Definition 2.3 (entropy process solution). An entropy process solution of (P) is 
a couple (/i, w) of measurable functions /i : Qt x (0, 1) ^ R and w : Qt R 
satisfying the following conditions: 

(DM) (regularity and consistency) pL £ L°°{Qt x (0, 1)), w E LP{0, T; W^J '^(f})), 

and (j){p{t, X, a)) = w{t, x) for a.e. (i, x, a) G Qt x (0, 1). 
(D'.2) For all ( G P([0, T) x fl), 

' b{p.)dt£, + \{p) ■ - a{n, Vw) • - i^ip)^ ] dxdtda 



+ / uo^iO,-)dx+ / / f^dxdt^O. 
Jn J Jqt 

(D'.3) For all (c, ^ G R=^ x V{[0,T) x H), ^ > 0, and also for all (c,^ G R x 
V{[0,T) X r!), ^>0, 

bf{p)dtC + qfiti) ■ Ve - (ritYMpyw) ■ Ve - iVcYil^Wl^)) dxdtda 

bt{uo)aO,-)dx+ f ff ivtYif^) f^ dxdtda >0. 
JoJ JQt 



Remark 2.3. In (D'.3), setting u := 

Jo 



fi{a) da one can rewrite the third term 



under the form 



(77±)'(M)a(Ai,V«;)= / / (r;±)'(u)a(ii, V^i;). 

Qt J JQt 

Indeed, we have w = 0(/i), and (f> is invertible on 'R\ E, so that fi(t,x,a) = 
(j>~^(w{t,x)) = u{t,x) whenever w{t,x) G R \ 4i{E); furthermore, Vw — a.e. on 
[w G 0(^5) U {0(c)} ], and the exact value of {rjfYip) on [ w G 4>{E) U {0(c)} ] does 
not matter because a{r, 0) = 0. For the same reasons, a(//, Vw) can be replaced by 
a(M, Vw) m (D'.2). 
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Remark 2.4. If u is an entropy solution of (P), then the couple {p-,w) defined 
a.e. on Qt x (0,1) (resp., on Qt) by ^{t,x,a) = uit^x) (resp., by w{t,x) = 
(j}{u{t,x))), is an entropy process solution of (P). In turn, if {lJ,,w) is an entropy 
process solution of (P) such that fi(t, x,a) = u{t,x) a.e. on Qt X (0,1) for some 
u : Qt ^ R; then the function u is an entropy solution of (P). 

2.3. Well-posedness of problem (P) in the framework of entropy solutions. 

First note the uniqueness result, which requires no range condition nor structure 
condition on the nonhnearities b and (j). 

Theorem 2.1. Assume that {IIi)-{IIr^) and {IIq)-{Hii) hold. 

(i) Assume that {fJ-,w) is an entropy process solution of [P). Then 

u{t,x)^ / ii{t,x,a)da 
Jo 

is an entropy solution of (P). Moreover, we have b{fi){t, x, a) = b{u){t,x) and 
tl)[pL)(t,x,a) = ip(u)(t,x) a.e. on Qt x (0,1). If {b + (f> + xp) is strictly increasing, 
then fi(t, X, a) = u(t, x) a.e. on Qt x (0, 1). 

(ii) Assume that u and u are two entropy solutions of (P) corresponding to the data 
Mo, / and iiQ, f , respectively. Then for a.e. t S (0, T), 

{b{u)-b{u))+{t)+ f ( {iiu)-^{u))+ 
n Jo Jn 

< f {b{uo) - b{uo))+ + I I sign+(M -«)(/-/). 
Jn Jo Jn 

(Hi) In particular, if u,u are two entropy solutions of (P), then b{u) = b{u) and 

ipiu) = ip{u). 

Remark 2.5. In Theorem 2.1(ii), on can replace u, resp. u, with an entropy 
subsolution, resp. with an entropy super-solution. The same proof applies. 

The following continuous dependence property is the central result of this paper. 

Theorem 2.2. Let {bn,'4'ny4>ni<^nifn',Uo, .fn), n G N, be a sequence converging to 
(6, ip, (j), a, f; uq, /) in the following sense: 

bm''Pm<t>n Converge pointwise to b,ip,(j) respectively; 
(4) • f„ , a„ converge to f , a, respectively, uniformly on compacts; 

^ri(uo) ^ ^{uq) in L^{n), and f„ —> f in L^(Qt). 

Assume that (6, -07 4>, a, f; mq, /) and {bn, ipn^ fpn, cin) fn; Uq, fn) (for eachn) satisfy the 
hypotheses (Hi), {H^), (H^), and {Hq)-{Hii), and, moreover, that the functions 
C(-) in (Hg), (Hio), and {Hn) as well as the i°°(0) and L^{0,T, L°^{n)) bounds 
in {H4) are independent of n. We denote by (P„) the analogue of problem (P) 
corresponding to the data and coefficients {bn, ipm (f'n, i^n, f Uq, fn). 

Assume that either &(R) — R, or the L°°[Qt) bounds on f^ in {II5) are inde- 
pendent of n. Assume that the structure condition {Hstr) holds, and (j) satisfies the 
technical hypotheses {H2),{II^). 

Let Un be an entropy solution of problem (P„). Then the functions ii„ converge 
to an entropy solution u of (P) in L°°{Qt) weakly-*, up to a subsequence. Further- 
more, the functions 4>n{un) converge to (^{u) in L^{Qt) up to a subsequence, and 
the whole sequences 6„(u„), ?/'n(w„) converge in L^{Qt) to b{u),'ip{u), respectively. 



WELL-POSEDNESS FOR DEGENERATE PARABOLIC EQUATIONS 



11 



Remark 2.6. The structure condition (Hstr) seems to be essential for results like 
Theorem 2.2 to hold (see [U]). 

The surjectivity assumption on b or assumption (H^) were only required to ensure 
the L°° estimate on u. One could work with unbounded solutions, in which case 
these assumptions can be replaced by a growth condition on f (see Remark 6.1 below). 

Notice that assuming simultaneously [Hstr), ^(R) = R- o,nd ijj = 0, by a change 
of the unknown u into v = b{u) we can always reduce the triply nonlinear problem 
(P) to the doubly nonlinear problem with b = Id. 

Finally, we state the well-posedness result for (P). Note that when only the data 
(mq,/) are perturbed, the continuous dependence result analogous to Theorem 2.2 
holds under the structure assumption {H'^^^) which is weaker than {Hstr)- 

Theorem 2.3. 

(i) Assume that {Hi)-{Hc,) and {He)-{Hii) hold. Then there exists an entropy 
solution to (P). Moreover, it is unique, in the sense of Theorem 2.1(iii). 

(ii) Assume in addition that the structure condition {H'^^^) holds. Then the entropy 
solution of (P) depends continuously on the data (uo,f). More precisely, let Un be 
an entropy solution of {P) with data {uQ,fn). Assume &(wq) — > 6(uo) L^{fl) and 
/„ — s- / m L^(Qt), as n oo. Assume that the bounds in {H4),{H^) are uniform 
in the sense 



f||L-(S2) < Const,- 



T 



• either b{+oo) = +00 and / \\f ^{t, dt < Const , or ip{+oo) — +oo 

Jq 

and ||/+||l-(Qt) < Const,- 

• either b{— oo) ^ —oo and / (i, •)||ioo(f2) < Const, or 7/'(—cx)) = —oo 

"'0 

and Wfnh^-iQr) < Const. 
Thenb{un),4'i'an) and (f){un) converge, respectively, tob{u),ip{u) and (j){u) in L^{Qt) 
as n oo, where u is an entropy solution of (P) with data {uq, /). 

Moreover, if we reinforce hypothesis (Hg) by requiring the uniform monotonicity 
of a(r, •) in the sense 

there exists C G C(R^;R'*') such that 
(a(r,0-a(r,7y)).(e-7y)>l/C(r,^), 

then also W(j>{un) converge to W(p{u) in {LP{Qt))^ and o(it„, V(/)(m„)) converge 
to a{u, \/(j){u)) in {L^ (Qt))'^ , as n oo. 

3. Notation and preliminary lemmas 

Let us give some notation which will be used throughout the remaining sections. 
We will use the notation like [u g P] for the sets like {{t, x) € Qt \ u{t, x) G P}. 
For a measurable set H, we denote by U// the characteristic function of H. For 
C R, we set 



Th{z) := / ^H{s)ds- 
Jo 



clearly, Th{-) is a Lipschitz function with Tff (0) = 0. 

We denote by G the image </'(P) by of the "exceptional set" E; recall that 
E is closed, G is closed and meas (G) — 0. We denote by J a generic open interval 
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in R \ and by J it image (/)(/) which is a generic open interval in R \ G. For all 
e > 0, we choose an open set Ge D G such that mcas (G^) < Const x e. We denote 
by the open set 0~-'^(Ge) which contains E. When {H^) holds, we can simply 
take Ge = G-^ := {z G R I dist (z, G) < e}. 

Now let us prove the representation property used in Definition 2.1. 

Lemma 3.1. Let ipe{-) be the function defined by 

ipg-.z^H^ I 9{s)dip{s), 
Jo 

for a continuous non- decreasing function if : R R and a bounded piecewise 
continuous function 6* : R R". Then there exists a Lipschitz continuous function 
(^9 : 'y9(R) — > R" such that for all z g R, 

Proof. If (p{z) — (p{z), then the measure dip{s) vanishes between z and z; thus 
ipeiz) — (pg{z) — / 9{s) dtf[s) is zero. Therefore Lpg is well defined. For all r,f £ 

J z 

(/j(R), (pe{r) - ipe{f) = Ve{z) - ^pe{z) = / 0{s)dip{s), where z e (p~'^{r), z e 

J z 

if^^{f). Thus 

\V9{r) - Mf)\ < \\e\\L^ \^{z) - ^{z)\ = \\9U^ \r-f\. □ 
Now, let us give a localized estimate of the gradient of w = 4>{u). 

Lemma 3.2. Let u be a bounded weak solution of (P). Then there exists a constant 
C depending on G(-) in (Hg), on ||w|U»(Q5,) and on ||fe(uo)||Li(J2), II/IIli(q) such 
that for all Borel measurable set _F C R, 



(5) J J[u<£F\ J J[w£<P{F)\ 

< CYsiTF(t>{-) = Gmeas(0(F)). 

Proof. Without loss of restriction, one can assume that F is bounded; indeed, 
otherwise we can replace F with F'^^ :— F O [— M, M] and then pass to the limit 
as M — !■ +00 in inequality (5) written for F^^ . 

Set H = (^(F), and note that Th{w) = ThWu)) € iP(0, T; '^(17)) can be 
approximated by admissible test functions in (2) of Definition 2.2; one has 

VTh{w) = VwUh(w) = V</>(u)llF(it), 

and 

(6) \\Th{-)\\.oo < f d(b{s)=VaTF(bi-)=meas{H). 

Jf 

Using this test function, with Remark 2.2 and the standard chain rule argument 
known as the Mignot-Bamberger and Alt-Luckhaus formula (see, e.g., Alt and 
Luckhaus [1], Otto [35], Carrillo and Wittbold [21]) we get 



(7) 



Bpiu){T, ■)+ ^{u)TH{<j>{u)) + / / a(ii, Vw;) • VT^H 

J JQt J JQt 

= ( Bpiuo) +11 fTniw) - / / fH • VT^(0(ii)), 

Jn J JQt J JQt 
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where Bp {z) = [ TH{(l){s))db{s). 
Jo 



The last term is zero thanks to the boundary 



condition w|s = 0. Indeed, because F is assumed bounded, f(-) is bounded on 
the support of T^{(p{-)); thus by Lemma 3.1 there exists a Lipschitz continuous 
vector- valued function such that 

f(s)drH(0(s)) = 0(</>(z)). 

Hence g o w e LP{0,T;WQ'^{fl;'R'^)), so that one can apply the Green-Gauss 
formula to get 

div ( / f(s) dTui^is))) ^ / / q{w) • = 0, 
1 Jo Jo Jdn 

where v is the exterior unit normal vector to dfl. By definition of Tff(-) and because 
(/>{■) is non-decreasing, dropping positive terms in the left-hand side of (7), by {Hg) 
we infer 

w) ■ Vw 



C(||u||l~(Qt)) J J[ueF] J J[w£ct,{F}] 

< {\\biuo)\\LHn) + \\f\\LHQ))\\TH{-)\\^. 
Hence the claim follows by (6). □ 

In the above proof, wc have used two chain rule lemmas. Now we notice that 
both apply for u{-) replaced with a "process function" a), as in (DM), provided 
that for a.e. a e (0, 1) one can substitute 

u{t, x) :~ / ii{t,x,a)da 



by fi{t, X, a) in the expression of the test function. 

Lemma 3.3. Let {^J.,w) satisfy (DM), and S : H ^ R be a Lipschitz continuous 
function such that S{0) = 0. Let C € L°°(0,T). Then 
fi 

f{n{t,x,a))) ■ WS{w{t,x)) C{t) dtdxda = 0. 



Proof. By Lemma 3.1, there exists a Lipschitz vector- valued function g such that 
f(s)d5(^(s)) = g(<?!)(z)) for \z\ < \\h\\l^{q^x{o,i))- We have 

'f(/.(a))da- V5(</)(u))C 



QtJO 



f(Ma))- V5(0(/x(a)))daC 
\iy { r^"\{s)dSWs)))daC 



[ [ div s{w) C = [ [ &iw) ■iyC = 
10 Jn Jo Jdn 



because for a.e. a e (0, 1), (t>{lJ-ia)) = 0(w) = w G LP{0, T; W^'^'ifl)). □ 
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Lemma 3.4. Let fl be a bounded domain 0/ R", T > 0, Qt ■= (0,T) x il, and 
1 < p < +00. Let g G C(R; R). Let b G C(R; R) be non- decreasing. Set 



Bg{z) := / g{s)db{s). 
Jo 

Let ^ e L°°{Qt X {0,1)); set u = / fi{a)da. As. 

Jo 



that 



g{u)eLPiO,T;W^^P{n))nL°^{QT) 

and, moreover, 

g{n{a)) = g{u). 

As.sume that 

dt(^J^ b{fi{a)) da) e LP'{0,T; W-''^' {n)) + L\Qt) 

and ^ 

/ b{fi{a)) da\t=o = b{uo) 
Jo 

in the following sense: 

e LP{0,T; Wo-P{n)) such that dt^ G L°°{Qt) and ^{T, ■) = 0, 



Proof (sketched). Note that the claim of Lemma 3.4 cannot be deduced directly 
from the usual Mignot-Bambcrgcr and Alt-Luckhaus chain rule lemma; the reason 
is that we cannot expect dtb{fi{a)) to belong to iP'(0, T; W~'^'P' (Q)) + L^{Qt) for 
a.e. a. But it suffices to reproduce the proof (see, e.g., [21]) which is by discretiza- 
tion of dt /q^ b{fi{a)) daj . Indeed, we have for a.e. t,t — h G (0, T), 
1 / .1 



{dt( 6(Ma))da) ,0 = - / / / b{fi{a))dadt^~ / 6K)C(0,-)- 
^Jo ' J JqtJo Jn 

Then for all ( e V{[0,T)), 

{dt ( / b{^l{a)) da) , 5(^)0 = - // / Bg{^^{a)) daCt~ f i?s("o)C(0). 



h 



b{fi{t, a)) da — / b{fi{t~h, a)) da \ g{u{t)) 



If 



= J - b{l^{t~h,a))) g{fi{t,a)) da, 

and now we can reason separately for each a. Thus the arguments of [21] apply. □ 

4. Proof of L^ contraction and comparison principles 

Now we turn to the proof of Theorem 2.1 and Remark 2.1. Most of the state- 
ments are standard. We only notice that while proving Theorem 2.1(i), one obtains 
that b{fi) and ipifJ-) &re independent of a; since 4'{fi) = w is independent of a by def- 
inition, one concludes that f(/.i) = ^(b{^),ip{^),4i{^)) is also independent of a, thus 
the entropy process solution gives rise to the entropy solution u = fJ-ia) da. 
This is the only point where the special structure of the dependency of f on m is 
used. 
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The proof of Theorem 2.1 is esscntiaUy the same as in CarriUo and Wittbold 
[21]; it is based on the techniques of CarriUo and on hypothesis (i?ii) (notice 
that in the case (j) = Id, one has E — 0; therefore (-ffn) reduces to the CarriUo- 
Wittbold hypothesis in this case). For Theorem 2.1(i), wc also need the adaptation 
of the CarriUo arguments to the framework of entropy process sohitions. This has 
been done by Eymard, Gallouct, Herbin and Michel [2-1], Michel and VovcUe [33] 
and Andreianov, Bendahmane and Karlsen [(i]. Therefore we only point out why 
hypothesis (Hn) is sufficient for the uniqueness of an entropy solution in the case 
of problem (P) with cf) that can be not strictly increasing. 

The role of hypothesis (i?ii) is to ensure that 

(8) limsup / / // i ( a(M, Vu;) - a(w, Vw) ) • ( Vu- - Vui) ll[o<^,_i,<e] > 0, 
where 

u= fi{a)da, w — (p{u), w = / jj,{a)da, w — 4>{u), 
Jo Jo 

and {n{t, X, a),w{t, x)) and (/i(s, y, a), 'w{t, x)) are two entropy process solutions of 

(P). Here, following Kruzhkov [28], we have taken two independent sets of the 

variables (t,x) and {s,y). 

We split the integration domain Qt x Qt into several pieces. 

First, notice that a.e. on [w e G] x Qt, we have Ww = 0; thus the integrand 
in (8) is reduced to a{u, Vw) Vw ]l[o<tii-ji)<E] : which is non-negative. The same 
argument applies on Qt x [w € G]. 

Thus it remains to investigate the integrand in (8) on the set [w ^ G] x [w ^ G]. 
Let us introduce G"^ := {z G R | dist (z, G) < e}. For a.e. {t, x, s, y) E [w ^ G]x[w ^ 
G] we have : 

(a) either 'w{t, x) and w{s, y) belong to the same connected component of R \ G; 

(b) or (a) fails, but w{t, x) eG''\G and w{s, y) £ G^ \ G; 

(c) or both (a) and (b) fail, but then \w{t,x) — w{s,y)\ > 2e. 

We then split [w ^ G]x[w ^ G] into disjoint union of sets 6*0 U 5^ U S'c, according 
to which of the above cases (a),(b),(c) takes place at {t, x, s, y) G [w ^ G]x [w ^ G]. 
On Sa, we use assumption (Hn) and infer that the integrand in (8) is minorated 

by 

-max{G(r,s) | |r|,|s| < ||u||oo} (1 + | Vw|p + | Vw|P) ll[o<„_A<e] • 

Because the 2(Af + l)-dimensional Lebesgue measure of the set [0 < w ~ w < e] 
goes to zero as e — > 0, the limit of the corresponding part of the integral in (8) is 
minorated by zero. 

On Sb, we bound the integrand in (8) from below by — | (| Vti;|'' + 1 Vw\p). Using 
Lemma 3.2 we have e.g. 

£ J J[weG^\Gy J[wi=G'\G] £ J J 

By the continuity of the Lebesgue measure and because U£>oG^ \ G = 0, the 
measure of the set [w e G'^\G] tends to zero as e — > 0. Therefore, using assumption 
{H3), we deduce that the corresponding part of the limit in (8) is non-negative. 

Finally, on Sc the integrand in (8) is zero. 

This ends the proof of (8). 
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5. A PRIORI ESTIMATES 

The following estimates are rather standard. 

Lemma 5.1. Let (6„, ipn, 4>m Am fni ' G N 6e a sequence of data satisfying 

the assumptions of Theorem 2.2. Assume that the limiting data (6, i/;, </>, a, f; uq, /) 
are such that (H^) and {Hstr) hold. 

Let Un he an entropy solution of problem {Pn). Then there exists a constant M 
and a modulus of continuity oj : R+ R+, such that for all n G N, 

(^) hnllL-(Q^) < M; 

(ii) the following quantities are all upper bounded by M : 

||(/>«(-U„)||^p(0^y.;^^i,P(j2)), ||(/>nKOIUp(QT)> \Wn{Un, V0„ (u„ ) ) II Lp' (Q^) , 
\\i^n{Un)(l)n{Un)\\L^(QT): 1 1 ^« ("n) 1 1 L~ (0,T;L1 (O)) 

where _B„ is defined in (1) with b,(j) replaced by 6„,0„; 

(Hi) for all A>0, // |</)„(M„(t+A, a;))-0„(u„(t, a;))| < uj{a). 

J JQj,_a 

Proof, (i) First assume 6(R) = R. Consider the function 

M(t) sup(||6„K^)||ioo(o) + / ||/«(r,-)||L-(n)rfT) < +00. 

nSN JO 

Then for any measurable choice oiu{t, x) € b~^{M{t)), u is an entropy supersolution 
of (Pn)- Similarly, u{t,x) G b~^{—M{t)) is an entropy subsolution of {Pn). The 
comparison principle of Remark 2.5 ensures that a.e. on Q, 

~M{T) < -M{t) < bniun){t,x) < M{t) < M{T). 

Now the assumption &(R) = R and the pointwise convergence of 6„ to h ensure the 
uniform L°°{Qt) bound on u„. 

If b{+oo) < +00, then {H5) ensures that any constant u g iW f^W l^^ (Qt)) 
is an entropy subsolution of (Pn)- As hereabove, the comparison principle and the 
pointwise convergence of ipn to "0 satisfying 0(+oo) = +00 yield a uniform majo- 
ration of u„. The case 6(— 00) > —00 is analogous. 

(ii) We use the test function (pniun) hi the weak formulation of {Pn)- The duality 
product between 

Mun) G LP{0,T;W^-P{n))nL'^{QT) 

and 

dtbn{un) G LP' {0,T;W-''P' {n)) + L\Qt) 

is treated via the standard chain rule argument ([1, 35, 21]). Using in addition the 
chain rule of Lemma 3.3, the L°° bound on m„ shown in (i), the uniform coercivity 
assumption (-ffg), and the obvious inequality Bn{z) < 6„(z)0„(z), we obtain the 
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inequality 

n ['^n j I 

n Jo Jn ^ 

< [ bn{u'f^) <l>n{Uo) + [ [ /„0„(U„) 

Jn Jo Jn 

with some c > independent of n. Notice that the functions 6„, 0„ are locaUy uni- 
formly bounded because they are monotone and converge pointwise to b, 4>, respec- 
tively. Therefore the right-hand side of the above inequality is bounded uniformly 
in n, thanks to (i) and to the uniform bounds on the data Uq and /„ in L°°(Q) and 
in L^{Qt), respectively. The uniform estimate of the left-hand side follows. We 
then estimate \\4'n{un)\\ lp(Qt) by the Poincarc inequality; the {L^ {Qt})^ bound 
on a„(w„, V0„(7i„)) follows from the growth assumption {Hiq). 

(iii) Let A > 0. The weak formulation of (P„) yields, for a.c. t,t + A ^ (0, T), 

ibn{Un)it + A) -5„(u„)(i))^ 

n 

t+A, 

{-fn{Un) + a{Un, 'V(f>n{Un)))-'V£, - '0„(u„)^ + 

n 

for all ^ S WQ^^ifl) n L°°{fl) (here and in the sequel, we drop the dependence on x 
in the notation). Wc take ^ = (f>n{un{t + a)) — and integrate in t, then 

use the Fubini theorem which makes appear the factor A in the right-hand side; 
with the estimates of (i),(ii) and the uniform bounds on the data, we deduce that 



(9) // |6„K)(< + a) -6„K)(t)| |0„K)(t-f A)-0„K)(t)| < d|A|. 

J JQt-A 

Here c? is a generic constant independent of n. In the sequel, denote by (r„)„ a 
generic sequence vanishing as n — s- cx). Notice that by the Dini theorem, 6„,0n 
converge to 6, 0, respectively, uniformly on compact subsets of R. By (i), it follows 
that we can replace 6„, (/)„ in (9) by 6, 4>, provided that a term r„ is added to the 
right-hand side of (9). 

Now notice that assumption (Hstr) implies that cf) := cj) o is a continuous 
function. Let M be the L°° bound for u„ in (i). Let tt be a concave modulus 
of continuity of <f) on [—A/, A/], 11 be its inverse, and n(r) = rll{r). Let tt be the 
inverse of H. Note that W is concave, continuous, 7f(0) = 0. Denote v = 6(w„)(i + A) 
and y = b{un){t). We have 

\^{v)~^{y)\^ [ I 7ffn(|0(«)-0(2/)|) 
Qt-a J JQt-a V 

< IOt-aI^^^t^ / / - 4>{y)\) 



\Qt-a\ 



T-A 



Since \cj){v) ~ <j){y)\ < Trdw — yl), we have W{\4'{v) — 0(y)|) < |w — y| and 

n(|,^(^) - 0(y)|) = nmv) - ~^{y)\mv) - 4>{y)\ < \v - y\ ^v) - ^(y)] 

= \b{Un){t + a) - b{Un){t)\ |</)K)(t + a) - (f>{Un)it)\. 
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Therefore the estimate (9) (with bn,4>n and dA replaced by 6, (/) and c/a 
respectively) implies 

|0(it„)(t + A)-<?!)(u„)(i)| 

T-A 

1 



< IQt-aItT — r / / \Un{t + A) - U„(t)| |0(u„)(i + A) - (/)(M„)(t)| 

|QT-A|7f( --^^ — tA) < d7f(dA + r„), 
\\Qt~\\ J 

and finally, replacing (j) with we get 



Now using the fact that r„ — > as n ^ oo and the fact that for all fixed n e N, 
the left-hand side of (10) tends to zero as A ^ 0, we deduce the claim of (iii). □ 

6. Proof of the general continuous dependence property 

In this section, we prove Theorem 2.2. First notice that the uniform estimates 
of Lemma 5.1 and Lemma 3.2 apply. It follows that there exists a (not relabelled) 
subsequence (u„)„ such that 

• Wn '■= (pniun) convcrgcs strongly in L^{Qt) and a.e. on Qt to some function w; 

• Wn converges weakly in LP{0,T;W^-^{fl)); 

• Xn ■= a(Mn, Vwn) convcrgcs weakly in {Qt) to some limit x; 

• Un converges to /i : Qt x (0, 1) — > R in the sense of the nonlinear weak-* 
convergence (3). 

Denote u(-) = /q ^{■,a)da. Thanks to the uniform L°° bound on m„ and to 
the uniform convergence of 0„ to on compact subsets of R, we can identify the 
limit of Wn{-) with <f>{ii{a, •)) da; moreover, 0(//(a, •)) is independent of a G (0, 1), 
because the convergence of 4>{un) to w is actually strong. Thus w, (j){u) and 4>{^{a)) 
coincide. We also identify the limit of Vw„ with Vui, because the two functions 
coincide as elements of V . 

The following lemma permits to deduce strong convergence of m„ to u on the set 
[u ^ E] (recall that E is assumed to be closed). 

Lemma 6.1. Let 4>n {-) be a sequence of continuous non- decreasing functions con- 
verging pointwise to a continuous function </>(•)• Assume (/){■) is strictly increasing 
on 'R\E. Let I be an open interval contained in R\E, and 0„(u„) — > a.e. on 
Q. Then u„ — > u a.e. on [u G /]. 

Proof. Let / — (a, b) and choose /' = (a', b') with a < a' < b' < b. Introduce d > 

by 

S := min{0(a') - 0(a), (j>{b) - (j){b')}. 
Notice that by the Dini theorem, the convergence of (f>n{') to (/){■) is uniform 
on all compact subset of R. Thus for all e > 0, there exists G N such that 
\\4'n — 011 c([a,h]) < £/2- With e ^ 6, it follows that for all ??. sufficiently large, 
(f>„{b) - (j){b') = (j)n{b) - (t)ib) + (f){b) - (i){b') > -6/2 + 5 ^ (5/2, and similarly, 
(j>{a') — 4>n{o-) > S/2. By the monotonicity of (/)„{■), 4>{-) , 

max |0„(z) - <j){u)\ = max{0„(6) - 0(&')i 't'ia') ~ </>«(«)} > (5/2. 
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Hence if u{t,x) G /' and if |0„(it„(t, x)) — (j){u{t,x))\ < 5/2, we have Un{t,x) G /. 
Since for all e > there exists N £ N such that x)) — (t){u{t,x))\ < e/2, 

we have in particular that for a.e. {t,x) G [u G /'], one has Un (z I for all n large 
enough. 

Thus for all e < S, for a.e. point of [u G /'] there exists N G N such that at this 
point one has 

|0(m„) - (j>{u)\ < \(f>{u„) ~ (t>n{Un)\ + - 

< Hn - 0||c([a,fc]) +£/2 < e. 

Therefore (fiiun) converges to 4'{u) a.e. on [u G /']. Since is continuously 
invertible on /, one also has Un u a.e. on [u G /']. Since / is open and /' d / is 
arbitrary, the claim of the lemma follows. □ 

Now we start to identify x with a{u, Vif). First, according to Remark 1.2, 
a{u, Vw) = a.e. on the set [u £ E] = [w & G]. Using Lemma 3.2, we now deduce 
that also x = on this set. 

Lemma 6.2. Let x be the weak (Qt) limit of the sequence 

Xn = an{Un, V(?!)„(u„)), 

and let Wn = 4>n{un) converge to w = (f'i'u) a.e.. Then x = a.e. on [u € E]. 
Moreover, Xn converges strongly to zero in L^([u £ E]). 

Proof. Since ll[tig£;] G L°°{Qt) C LP{Qt), by the definition of the weak conver- 
gence, the function x^ueE] = X^weG] is the weak {Lp' (Qt))^ limit of Xnll[toGG]- 
For all e > 0, choose an open set G^ D G oi measure less than e. Because Wn ^ w 
a.e. and G^ is an open neighbourhood of G, we have 

[u e E] = [w e G] = RU (Ujvgn[«^ G G, w„ G GeVn > N]), 
where meas (i?) = 0. Since 



[w G G, Wn G GeVn > N] 

■ NeN 

is an increasing collection of sets, the corresponding sequence of measures converges 
to meas ([w G G]) as ^ oo, by the continuity of the Lebesgue measure. Because 

meas ([if G G]) > meas {[w G G, wn G G^]) > meas {[w £ G, Wn £ G^ V?i > iV]), 

we conclude that meas([w £ G, Wn ^ G^]) tends to zero as n ^ oo. Now by the 
Holder inequality. 



Ix«11[«,gg]I|li(Q) = / / |a„(u„, Vw„)| 

J J[wGG] 
< / / \an{Un, Vw„)| 



Qt 



meas {[we G, Wn ^ G^]) 

For all fixed e > 0, the last term tends to zero as ri — > oo, thanks to the boundcdness 
of the sequence a„(u„, Vw„) in {Qt)- The first term in the right-hand side 
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converges to zero as e uniformly in n; indeed, by (-ffg) and by Lemma 3.2, it 
is majorated by -!-mcas(Ge) < s/c 

It follows that = a.c. on Qt, which ends the proof. □ 

Now we use the Minty-Browder argument to identify x with a(u, Vw) on the 
sets [m S /], for all open / C R \ -E. A crucial role is played by Lemma 3.3, which 
permits us to pass to the limit in the product 

of two weakly converging sequences. 

We proceed in the classical way, but use the test functions Tj{w) (at the limit) 
and Tj{wn) (before the passage to the Hmit), where J — and 

Tj{z) := [ Aj{s)ds 

JO 

is the truncation function that localizes the values of the solution to the inter- 
val /. Notice that we can assume that VTj(u'„) converges to VTj(w) weakly 
in (LP(Qt))^. Indeed, the sequence (Tj(w„))^ is bounded in LP{0,T;W^'''{n)), 
hence (up to a subsequence) it converges to a limit; this limit is identified with 
Tj(w) since w„ converges to w in L^{Qt), and Tj is Lipschitz continuous. Now let 
us give the details. 

We first pass to the limit into the weak formulation of (P„) (see Remark 2.2), 
getting 



(11) 



(9t^y b{^{a)) da^ + div^ J f{^{a))da^+ J ^{^(a)) da ~ div x + / 

mLP'iO,T;W-^^p'in))+L\QT), 



1 

b{^{a)) da\t=a = b{uo) 



in the same sense as in Remark 2.2. Take ^ G X'([0,T)). Because w 
(f>{fj,{a)) for a.e. a € (0, 1), using Lemma 3.4 we have 

j^(dt[j\[ti{a))da) , T,{<P{u))(: 

/ BJ{^l{a))da^t(:- I B.j{uo)m 
iqtjo Jn 

with Bj{z) / Tj{<j>{s)) db{s). Define similarly S7(z) := / Tj{(pn{s)) dbn{s). 

Jo jo 
By the standard chain rule lemma of [1, 35, 21] we get 

(13) r{dMUn),Tj{MUn))0 = - I I B"j{Ur,)dt(:~ I B^jiu^)^. 



One shows easily that i?" converges to Bj uniformly on compact subsets of R, 
because of (4). In particular, Bj{uq) converge to Bj{uq) in L^{U). Moreover, the 
nonlinear weak-Tk- convergence of Un to n yields 



lim // B^j{u^)dt(:^ Bj{ti{a))dadtC. 

Thus the left-hand side of (12) coincides with the "n +oo limit" of the left-hand 
side of (13). 
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Similarly, the nonlinear weak-* convergence of u„ to fi permits to conclude that 
da Tj{4>{u))C ^ f f da ( 

= lim / / i!n{Un)Tj{(j)n{Un))C- 

Now let us take the test functions Tj{wn)C ^-nd Tj{w)C in the weak formulation 
of (P„) and in (11), respectively. Without loss of restriction, we can assume that 

t = r is a Lebesgue point of / Bj{ii{a)) da. Using the above calculations and 

Jo 

Lemma 3.3, then letting ({t) converge to 11[o,t)(^) we deduce the equality 
(14) // X - ^Tj{w) - lim / / a„(ii„, V«;„) • WTj{wn). 



n — >oo 



Lemma 6.3. Assume that for all E R^, a,i(-,^) converge to o(-,^) uniformly on 
compact subsets of an open interval I C R\ With the notation and assumptions 
above, one has x = fl('", Vw) a.e. on [u G /] = G J]. 

Proof. Take an arbitrary function C G {L'p{Qt))^ such that C = a.e. on [w ^ J]. 
Take A G R. Set C = ^Tj{w) + X(; we have C = a.e. on [w ^ J]. By the classical 
Minty-Browder argument, considering ±A j one concludes that x = Ww) 
a.e. on [w G J], provided the following relations can be justified: 



X ■ {VTj{w) - C) > hminf / / a„(M„, Vu;„) • iVTj{wn) - C) 

J Jqt 

> liminf / / o„(m„,C) • (VTj(u;„) -C) 



n^oo 



(15) 

= liminf// a(u,C)-(VTj(u.„)-C) 



a{u,C)-{VTj{w)-C). 

Now we justify (15). Because (14) holds and WTj{wn) converges to WTj{w) weakly 
in {LP{Qt))'^ , the first inequality and the last equality in (15) are clear. 

The second inequality comes from the monotonicity of a„(M„, •). Indeed, by the 
choice of 

o„(u„, Vw„) • ( VTj(w„) - C) 
a„(M„, Vw„) • Vw„ 



>0 = / / a„(?/„,C)-(Vrj(u;„)-C), 

because a„(-,0) = 0. Further, since J is open and u)„ w a.e. on Qt, as in 
the proof of Lemma 6.2 we have meas([w G J,Wn ^ J]) ^ as n — > oo. On 
[w G J,Wn G J] one has a„(u„, Vw„) = a„(u„, VTj{wn)), and (ifg) applies. On 
[lu G J, tfn ^ J] the terms with VTj{wn) are zero; finally, as in the proof of 
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Lemma 6.2, the integrals 



a„(M„, Vw„) • C and / / a„(M„,C) • C 

[w£j,w„^J] J J[w£j,w„^J] 

both tend to zero, by the equi-integrabiHty argument. 

It remains to justify the last but one equality in (15). Thanks to Lemma 6.1, we 
have Un — > w a.e. on [w G J] = [u G I]. In particular, a.e. on [u £ /] one has Un & I 
for sufficiently large n. Using the locally uniform on / convergence of a„(-,^) to 
o(-,^), we conclude that a„(w„,C) converges to a(u, Q a.e. on [u € /]. By (i/io) 
and the Lebesguc dominated convergence theorem, (15) follows. □ 

Finally, notice that thanks to (DM) and Remark 1.2, we have a(/x, Vw) = 
a{u, Vw). The identification of x thus being complete, from (11) we readily con- 
clude that (fj.jw) verifies (D'.2). 

Now we can pass to the limit in the entropy inequalities corresponding to (Pn) 
and deduce (D'.3). Because regular boundary entropy pairs {ilt.sT^t.e) "^^n be ap- 
proximated by convex combinations of scmi-Kruzhkov pairs, we have the analogue 
of (D.3) for Un with (?7^, q^) replaced by (?7^£, i^c^^) (with q"'^''^ converging to q^ 
uniformly on compact subsets of R) . 

Consider the third term in (D.3). We have 



lim // (ry^^)'(u„)a„(u„, V-u;„) • 

(16) ""^^JJqt ' 

= lim ( / / + / / )(?y±J'(w„)V^a„(u„,V«;„). 



ueG] J J[w<^G\/ 

By Lemma 6.2 and because {rf^e)' {un)"^ S, are bounded uniformly in n, the first 
term converges to zero; also notice that 

0= // (7y±J'(7.)Ve.a(w,VH 

J J[weG] 

because a(u,Vw) = a(it,0) = a.e. on [w e G\. By Lemma 6.1, we have w„ — > u 
a.e. on [w ^ G]; by the dominated convergence theorem, we infer that {T]f^)'{un) 
converges to (?7^e)'(u) strongly in L^{[w ^ G\). Because also Xn = a,n{un,^Wn) 
converges to x = a(u,Vu;) weakly in LP {[w ^ G]), the second term in the right- 
hand side of (16) converges to 

(^c*.)'(^)Ve-a(u,Vu;). 

w^G] 

The passage to the limit as e ^ in the other terms in (D.3) is straightforward, 
using the uniform boundedncss of {un)n and the nonlinear weak-* convergence 
property (3). At the limit, we conclude that (D'.3) also holds. Thus (/i, w) is an 
entropy process solution of (P). 

Now by the result of Theorem 2.1(1), {fi,w) gives rise to an entropy solution 

u ~ m(q^) '^Q^ 
Jo 

of (P). By (3) (with F = Id), we conclude that (m„)„ possesses a subsequence that 
converges in L°°(Qt) weak-Tk- to u; we have already shown that the corresponding 
subsequence (j)n{un) converges to 4>{u) in L^{Qt)- 
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Moreover, b{fi{a)) and ?/'(/i(a)) are in fact independent of a. By Theorem 2.1(iii), 
we also have the uniqueness of b{u) and tp{u) such that u is an entropy solution 
of (P). By the well-known result of the nonlinear weak-* convergence (see e.g., 
Hungerbiihler [ -' ]), it follows that the whole sequences (bn{un))n and (jpn{un))n 
converge to b{u) and respectively, in measure on Qt and in L^{Qt)- 

This ends the proof of Theorem 2.2. 

Remark 6.1. In the case assumption (H^) is dropped, in order to deduce that u is 
an entropy process solution, along with the assumption (Hstr) one needs a growth 
restriction on f of the following kind: there exists a function M £ C(R+; R^) and 
a function C with lim C{z)/z = such that 

z — ^ + oo 

|f(6W, <A(r), ^(r))| < M(|6(r)|) £ (^\c^{r)\P + cp{s) db{s) + V'W^W 

for all r £ R, and the same assumption with |f(6(r), 0(r), 'ip{r))\ replaced with \b{r)\ 
and with \'ip{r)\. Indeed, these inequalities make it possible to use the nonlinear 
weak-k convergence framework of Ball [8] and Hungerbiihler [!()] without the uniform 
L°° bound on {un)n- 

7. Proof of Theorem 2.3 

In this section, we prove Theorem 2.3. The uniqueness claim was shown in 
Theorem 2.1; also notice that the continuous dependence result under the structure 
assumption (Hstr) was shown in Theorem 2.2. Let us first prove the existence claim. 

(i) First, consider the case where assumption {Hgtr) is fulfilled. Consider an 
approximation of (P) by regular problems with data ip, a, f; mq, /) such 
that the assumptions of Theorem 2.2 are fulfilled, and 6„, [&n]~"'^, f™, </>«, ['/>n]~^ are 
Lipschitz continuous on R. Using classical methods (cf. Alt and Luckhaus [I], 
Lions [Hii]), one shows that there exists a weak solution it„ G L'p{Q,T]Wq'^{VI)) to 
the problem {Pn) in the sense 

dtbn{u„) + div f„(u„) + V'(u,i) = div a{u„, V0„(u„)) -I- / 

in LP'{0,T; W-^^p' (n)) + L\Qt), with initial data 

bniUn)\t=0 = bn{uo). 

In addition, u„ verifies the entropy formulation (D.3), obtained along the lines of 
Carrillo [20]. By Theorem 2.2, we conclude that there exists an entropy solution of 

To prove existence without the structure condition (Hstr), we use the particular 
multi-step approximation approach of Ammar and Wittbold [4] (see also Ammar 
and Redwane [•■)]). We replace b by bk '■= ^+^Id and ij}, by V'm,n := V'+^Id^ + ;^Id^. 
The result of (i) for the corresponding problem (P^ is already proved. 

There exists a function ujjj „, constructed by means of the nonlinear semigroup 
theory (see e.g., [12]), such that 6/c(wJ^„) is the unique integral solution to the 
abstract evolution problem associated with (P^ (here and below, we refer to 
Ammar and Wittbold [4], Ammar and Redwane [3] for details). One then shows 
that „ coincides with the unique entropy solution of (P,^ „), the existence of this 
entropy solution being already shown. Further, the whole set [u^ ^)k,m,n verifies 
the uniform a priori estimates of Lemmas 5.1, 3.2. 



24 



B. ANDREIANOV, M. BENDAHMANE, K. H. KARLSEN, AND S. OUARO 



We then pass to the hmit in „ in the foUowing order: first k +00, then 
n — > +00, m +00. 

While letting k +00, we use the fact that tpm^ Lipschitz continuous. The 
fundamental estimates for the semigroup solutions permit to show that tpm,nium.n) 
are uniformly continuous on (0,r) with values in L^{fl); thus we get the strong 
precompactness of (uj^ „)fc in L^{Qt)- Thus, up to a subsequence, uj^ „ converge 
to Um,n which is an entropy solution of problem {Pm,n) corresponding to the data 

f; uo,f). 

Finally, we use the inequalities Um+i,n < Um,n l£ Um,n+i which follow readily 
form the comparison principle of Theorem 2.1(ii). The monotonicity argument 
yields the strong convergence of Um,n- Now the whole scheme of the proof of 
Theorem 2.2 applies with considerable simplifications, because no nonlinear weak-* 
convergence arguments are not needed. Passing to the limit in Um,n, we conclude 
that the limit u is an entropy solution of the original problem (P). This ends the 
existence proof. 

(ii) Existence for the limit data (uq,/) is now shown and we can apply Theo- 
rem 2.1 (ii). Then we deduce the L^((3t) convergence of 6(m„), '0(u„) to 
respectively. The convergence of 0(m„) to (f>{u) in L^{Qt) follows, by hypothesis 
{H'g^^) and because our assumptions imply the uniform L°°{Qt) bound on m„. 

The remaining claim of the strong (L^ (Qt))^ convergence of V4){un) is a rather 
standard part of the Minty-Browder trick. The argument is based upon the proof 
of Theorem 2.2. Because we already have the strong compactness of ((/)(«„))„, we 
can bypass the hypothesis {Hgtr) and deduce that the L°° weak-* limit u of u„ is 
an entropy solution of (P) with the data (ito,/). In particular, the {LP (Qt))^ 
weak hmit x of a(u„, V(f){un)) is equal to a(u, V(^(m)). Now notice that by The- 
orem 2.1(iii), under the structure condition (H'^f.^.) we also have the uniqueness of 
4>{u) such that u is an entropy solution of (P); moreover, by Remark 1.2, we also 
have the uniqueness of a{u, V0(u)) such that u is an entropy solution of (P). Thus 
X coincides with a{u, V(^(m)), so that (14) now reads as 

(17) // a(u,V0(M)) • V</)(u) = lim // a„(u„, Vw„) • V0(u„). 
It follows by the weak convergences of W(j){un) and of a(w„, V4>{un)) that 

(18) lim / / {a{u, V0(u)) - a(u„, V0(u„))) • ( V0(u) - V0(u„)) - 0. 

Notice that a.e. on the set [u G P], V(/)(u„) converges to = W<f>{u), by Lemma 6.2 
and by the coercivity assumption (Pg). On the set [u ^ P], we can use Lemma 6.1 
to replace a{u,\/cj){u)) with a(u„, V(/)(u)) in the above formula (18). Using the 
uniform monotonicity assumption (Pg), we can now conclude that the convergence 
of V(/)(m„) to V0(u) holds a.e. on Qt- 

Separating again the sets [u e E] and [u ^ P], we deduce that the sequence 
of nonnegative functions a(w„, V(/)(u,i)) • \7<f){un) converges to a(w, V(/)(u)) • \i'4>{u) 
a.e. on Qt- Together with (17), this implies that a(it„, V0(u„)) • V4){un) also con- 
verge to a(u, V(/)(u)) • V0(ti) in L^{Qt)', in particular, they are equi-integrable on 
Qt- The coercivity assumption (Pg) now implies the equi-integrability on Qt of 
the functions | V(/)(m„)|p. Combining this argument with the a.e. convergence of 
V(/>(Mn), we deduce our claim from the Vitali theorem. The (P^ {Qt))^ conver- 
gence of a{un, S/cpiun)) to o(u, \7(j){u)) follows in the same way. 
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8. WELL-POSEDNESS for the DOUBLY NONLINEAR ELLIPTIC PROBLEM 

Wc first notice tliat tlic wcll-poscdncss result for tlic degenerate elliptic problem 



[ M = on dil 
(i?^) s € and V'(R) = R 

follows from Theorem 2.3, upon setting & = 0, f{t, •) = s(-) and arbitrarily prescrib- 
ing uq. The definition of an entropy solution of (S) can also be formally obtained 
from Definition 2.2. 

Let us notice that the analogue of the general continuous dependence property 
of Theorem 2.2 holds without any additional structure condition: 

Theorem 8.1. Let (ipm 4>m cin, fn; s), n Cz be a sequence converging to (ijj, (f>, a, f; s) 
in the following sense: 



Assume that ip, a, f; s) and {ipn, ^nj c„, f„; s„) (for each n) satisfy the hypotheses 
{Hi), (Hq)-(Hii), and (H^), and, moreover, that the functions C in (Hg), (Hiq), 
and (Hii) as well as the L°°(D,) bound in (H^) are independent of n. We denote 
by (Sn) the analogue of problem (S) corresponding to the data and coefficients 



Assume that (p satisfies the technical hypotheses (H2) ,{H^) . 

LetUn be an entropy solution of problem (Sn). Then the functions u„ converge to 
an entropy solution u of {S) in L°°(D,) weakly-*, up to a subsequence. Furthermore, 
the functions (j)n{un) converge to 4>{u) in LP(il) up to a subsequence, and the whole 
sequence ipniun) converges to ^^{u) in L^(^l). 

The proof of Theorem 8.1 is contained in the one of Theorem 2.2, because the ^^(ri) 
bound on Vwn is sufficient for the strong precompactness of Wn- 




ipmfpn converge pointwise to ■0,0 respectively; 

f„,a„ converge to f , a, respectively, uniformly on compacts; 

Sn converges to s in L^{Q,). 



(V'ri; 071; ^ri: f n 5 ^: 
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